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Abstract 

This paper provides results for the spectra of triply charmed and bottom baryons based on a 
constituent quark model approach. We take advantage of the assumption that potential models 
are expected to describe triply heavy baryons to a similar degree of accuracy as the successful 
results obtained in the charmonium and bottomonium sectors. The high precision calculation of 
the ground state and positive and negative parity excited states recently reported by nonpertur¬ 
bative lattice QCD provides us with a unique opportunity to confront model predictions with 
data. This comparison may also help to build a bridge between two difficult to reconcile lattice 
QCD results, namely, the lattice SU(3) QCD static three-quark potential and the recent results of 
nonperturbative lattice QCD for the triply heavy-baryon spectra. 

TAGS numbers: 14.40.Lb,12.39.Pn,12.40.-y 
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I. INTRODUCTION 


In a previous publication [l| we have pointed out that the static three-quark potential 
with parameters determined from SU(3) lattice QCD does not reproduce the triply heavy- 
baryon bbb and ccc spectra measured also in lattice QCD We argued several possible 

reasons for such disagreement. In this work we aim to analyze whether the triply heavy- 
baryon spectra recently calculated by means of nonperturbative lattice QCD techniques can 
be understood within a constituent quark-model framework by means of simple Cornell-like 
potentials as it is the case of charmonium and bottomonium spectra This calculation 

may help to connect the description of the heavy-baryon spectra obtained by nonperturbative 
lattice QCD techniques and the static potentials derived within SU(3) lattice QCD. 

Theoretically, one would expect that potential models would be able to describe triply 
heavy baryons to a similar degree of precision as their success in charmonium and bottomo¬ 
nium. As noticed by Bjorken some time ago j^, bound states of three-heavy quarks reveal a 
pure baryonic spectrum without light-quark complications. In the same way the QQ inter¬ 
actions are examined in heavy mesons, the study of triply heavy baryons will probe the QQ 
interactions in the heavy quark sector. However, no experimental results are available so far 
for triply heavy baryons (see Ref. [l^ for a recent calculation of production cross section 
at the LHC), and thus, the predictions of their properties cannot yet be compared to the 
real world. The recent precise calculation of the ground and excited states of triply bottom 
baryons [1, , together with the ground and excited states of triply charm baryons i , pro¬ 

vides us with a unique opportunity to test phenomenological quark models for baryons in 
the energy regime in which the description using potential models is expected to work best, 
the heavy-quark sector. Hence, the quark-model dependent calculations could be tested by 
comparing them to nonperturbative hrst-principles calculations in lattice QCD of the bbb 
and ccc systems. 

The paper is organized as follows. In the next section we will briefly described the Cornell 
potential used for the description of the heavy-meson sector. We will try to connect the well 
established parameters used in the literature for the heavy-meson sector into the relatively 
unknown heavy-baryon sector. In Sec. IHIl we will present our results. We will hrstly analyze 
the pattern of the nonperturbative lattice QCD results looking for the general structure of the 
potential. We will pursue different parametrizations in a trial to get a unihed description 
of the bbb and ccc spectra. The results will be derived by different numerical techniques 
previously tested by our group: generalized Gaussians variational approaches, hyperspherical 
harmonics and Faddeev equations [IM3. Finally, in Sec. IIVI we will summarize the main 
conclusions of this study. 


II. A CONSTITUENT QUARK MODEL POTENTIAL FOR THREE-HEAVY 
QUARKS. 


Since the early days of QCD the interaction among heavy quarks has been explored as 
an important tool to learn about the behavior of QCD at low energies. At the end of 1974, 
when the new particles seen at Brookhaven and SLAC were identihed as cc bound states, 
explicit models were proposed to calculate the spectrum and the radiative transitions (see 
Ref. [l^ for a recent pedagogical review). The potential proposed in 15| is known as the 
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funnel or Cornell potential and it reads, 

Vg<i(r) = “ 


n - Vi 


+ b \fi — fA + c. 


( 1 ) 


Solving the two-body problem, one can tune the parameters to reproduce the low levels of 
charmonium. This was done by several groups in the 70’s and the authors were able to 
predict the missing states. The game became more challenging when the first bottomonium 
levels were found, trying to reproduce simultaneously the cc and bb spectra. Indeed, even so 
the interquark potential was not derived from QCD in early quarkonium phenomenology, it 
was assumed to be universal, or flavor independent. In QCD the gluons couple to the color, 
hence it is reasonable to assume that the potential is flavor independent. 

The QQ potential has also been extensively studied by lattice gauge theories [l^, being 
nowadays a very well-known quantity resembling the structure derived from the heavy-meson 
spectra. Thus, the typical shape of the color-singlet QQ potential is characterized by a short- 
range Coulomb behavior and a long-range linear rise, that well represents the double nature 
of QCD as an asymptotically free and infrared confined theory. The excitation spectrum 
of the gluon field around a static quark-antiquark pair has also been explored by lattice 
calculations [T^. On the large length scale the spectrum agrees with that expected for 
string-like excitations while in the short range it shows a Coulomb-like behavior as it was 
hrst noted within the context of the static bag picture of gluon excitations fl^ . 

The 3Q potential should be the analog of the famous Cornell potential for quarkonium. 
The short-distance behavior of V^^{r) is expected to be described by the two-body Coulomb 
potential as the one-gluon exchange result in perturbative QCD. It should be extended for 
the baryon case, with a factor 1/2 in front of its strength due to color factors [l^. As for the 
QQ case, the characteristic signature of the long-range non-Abelian dynamics is believed to 
be a linear rising of the static interaction. Moreover, the general expectation for the baryonic 
case is that, at least classically, the strings meet at the so-called Fermat (or Torricelli) point, 
which has minimum distance from the three sources (D-shape configuration) (To], 2^. The 
conhning short-range 3Q potential could be also described as the sum of two-body potentials 
(A—shape or linear configuration) [l9-22|. We have demonstrated in Ref. [l| the equivalence 
of both prescriptions for the case of triply heavy baryons (see Table II of that reference) and 
for different values of the heavy-quark mass. Thus, a minimal model to study 3Q systems 
may come given by. 






+ c. 


( 2 ) 
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The value of the QQ confinement strength b is usually fixed to reproduce that obtained 
from the linear Regge trajectories of the pseudoscalar vr and K mesons, = (429±2) 


MeV [16|. In the case of baryons, the linear string tension B is considered to be of the order 


of a factor 1/2 of the QQ case. The reduction factor in the string tension can be naturally 
understood as a geometrical factor rather than a color factor, due to the ratio between 
the minimal distance joining three heavy quarks and the perimeter length of a 3Q triangle, 
suggesting B = (0.50 ~ 0.58) b j^. For the particular case of quarks in an equilateral triangle 
B = ^b = 0.58 6 [i^. When the linear ansatz is adopted for the two-body potential, still 

the same relation holds for the strength of the Coulomb potential A ~ |a, due to color 
factors. The A ansatz (linear potential) has been widely adopted in the nonrelativistic 
quark model because of its simplicity jH, 13. 23-^ 
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III. RESULTS AND DISCUSSION 


To check whether the Cornell-like potential of Eq. ([2]) reproduces the bbb and ccc baryon 
spectra measured in lattice QCD we will make use of three different numerical meth¬ 
ods: the generalized Gaussians variational approach [ll|, hyperspherical harmonics jl^ and 
Faddeev equations [l3[. The three methods have been used and the difference in results is 
negligible. In all cases we solve the nonrelativistic Schrodinger equation 

{Ho + V^^{r)} T(f) = ET(f), 

where Hq is the free part of three-heavy quarks without center-of-mass-motion 

and Mq is the mass of the heavy quark. The mass of the heavy baryon will be finally given 
by Moq = ^Mq + E. 

It has been demonstrated in Ref. [l| that the results of nonperturbative lattice QCD for 
the bbb and ccc systems are not reproduced by the static three-quark potential derived in 
SU(3) lattice QCD j^. We had hypothesized on several possible reasons for this disagree¬ 
ment. Among them, we suggested the possibility that the parameters obtained in quenched 
QCD in Ref. might be different from those in 2-1-1 flavor QCD employed in Refs. 

To check this possibility parameters in 2-1-1 flavor QCD has to be obtained, unfortunately 
such calculations do not exist so far. We also wondered about the possibility that the fitting 
form Eq. (j2]) might not be appropriate to describe the static three-quark potential in lattice 
QCD or that the quark model description with ” 3-quark potential” might be inappropriate 
for bbb and ccc systems. We will try to mark out these possibilities by studying the nonper¬ 
turbative lattice QCD data in terms of the suggested Cornell-like potential of Eq. ([2]) with 
parameters adjusted to the data, looking also for the possible need of higher order terms. 

In Ref. [I| it has been noticed how bbb nonperturbative lattice QCD calculations point 
to a 3Q potential given by a mixture of a linear confinement and a Coulomb interaction. 
The comparison performed of the bbb and ccc nonperturbative lattice QCD results and 
calculations based on a mixture of a A-shape confinement and a Coulomb interaction using 
the parameters reported in Ref. i [A = 0.1410 and B = 0.0925 GeV^) showed a large 
difference in the excited states, predicting a small splitting between positive and negative 
parity excited states and also a small excitation energy for the positive parity states. These 
results may point to a lack of strength either in the confining or in the Coulomb potential. 

We have analyzed which of the two terms, linear confinement or Coulomb, would play the 
relevant role by keeping the strength of one of them as predicted by SU(3) lattice QCD and 
increasing the other to pursue a spectrum close to the nonperturbative lattice QCD results 
of Refs, [sl-lst]. We show in Fig. [1] results for the bbb system. In the upper panel we have kept 
constant the strength of the linear confining interaction, B = 0.0925 GeV^, and we have 
augmented the strength of the Coulomb potential up to A = 0.2787, a value rather close 
to the QQ Coulomb strength j^. Alternatively, in the lower panel we have kept constant 
the strength of the Coulomb potential A = 0.1410, and we have increased the strength of 
the linear confining interaction up to R = 0.1517 GeV^, which is also rather close to the 
QQ strength [^. One can easily identify the peculiarities of three identical heavy-quark 
baryons. By increasing the Goulomb potential, upper panel, the mass of all excitations is 
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FIG. 1: bbb excited state spectra, solid lines, for the potential of Eq. ([2]) with A = 0.2787 and B = 
0.0925 GeV^ (upper panel) or B = 0.1517 GeV^ and A = 0.1410 (lowe^anel). See text for details. 
The boxes stand for the nonperturbative lattice QGD results of Ref. [^. 


increased but the mass difference between the positive and negative parity excited states, 
E2(3/2~^) — Ei(l/2~), is smaller than predicted by nonperturbative lattice QCD. This mass 
difference could be enlarged by increasing the strength of the linear conhning potential as 
we have done in the lower panel, but without lowering too much the mass of the negative 
parity states, Ei(l/2~) and Ei(3/2~). It seems therefore that a moderate combined increase 
of the strength of the two terms. Coulomb and conhnement, may allow to get a better £t to 
the nonperturbative lattice QCD data. 

Thus, we have found a nice £t of the nonperturbative lattice QCD results with the poten- 


5 
















JP 



J'’ 

FIG. 2: bbb (upper panel) and ccc (lower panel) excited state spectra, solid lines, for the potential 
of Eq. ([2]) with the parameters A = 0.1875 and B = 0.1374 GeV^. The boxes stand for the 
nonperturbative lattice QCD results of Ref. for the bbb system and Refs. [^, 281 for the ccc 
system. 


tial of Eq. ([2]) and the parameters A = 0.1875 and B = 0.1374 GeV^. The resnlts are shown 
in Fig. |2]for the bbb and ccc systems. At this point one shonld bear in mind that the descrip¬ 
tion of the charmoninm and bottomoninm experimental data with a Cornell-like potential 
needs of a large valne of the QQ Conlomb strength, a ~ 0.51 — 0.52 [ol-lsj, as compared to 
that obtained in SU(3) lattice QCD j^. Using the phenomenological valnes reprodncing the 
experimental data of the heavy qnark meson spectra and the nonpertnrbative lattice QCD 
resnlts for the triply heavy baryon spectra, one wonld conclnde that A/a < 1/2, slightly 
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TABLE I: Excited states, in MeV, of ccc and bbb systems from nonperturbative lattice QCD, 
Refs, js, 2S| and Ref. and results for the ALl potential model of Ref. 25|. 



Ref. [2^ 

ccc 

Refs. [5,^ 

Ref. [2^ 

bbb 

Ref. 

Mf) 

-EiiV) 

472 

554±31 

451 

469±14 



567 

667±13 

581 

587±18 


-Eiil^) 

542 

636±13 

555 

567±18 

Elif) 

-Eiil^) 

316 

357±9 

360 

335.5±9.5 


different from 1/2 as the one-gluon exchange result. This conclusion could have been already 
anticipated from models designed to describe the light and strange baryon spectra and later 
extrapolated to triply-heavy baryons. The ALl model of Ref. [2^ made use oi A = 0.2534 
and B = 0.0827 GeV^. We show in Table [H the triply heavy baryon spectra obtained with 
this potential model as compared to the recent nonperturbative lattice QCD results. We 
observe how for the bbb system, it drives to a too high negative parity states together with a 
small splitting between positive and negative parity states (a similar situation to the upper 
panel of Fig. [1]) as a consequence of the large strength of the Coulomb potential. Besides, 
for the ccc system this model drives to rather low excited states as a consequence of the 
small strength of conhnement. 

The description of the bbb system shown in Fig. [2] has been obtained with a reasonable 
mass for the b quark, m;, = 4.655 GeV, however in the case of the ccc system one needs 

a large unrealistic mass for the c quark, rric = 2.050 GeV, what would drive to a ccc 

ground state with a mass around 7 GeV. This seems to indicate that triply charm baryons 

present a more involved structure that does not fit so nicely in a simple pairwise linear 

plus Coulomb interaction. Such deviation may has partially its origin on the fact that the 
lattice QCD calculation of Ref. did not address all the systematic uncertainties, but 
only statistical uncertainties are given. In fact, a calculation of the charmonium spectrum 


with the same lattice action and the same lattice spacing can be found in Ref. [2^ where 


one can get an idea of the typical size of the systematic uncertainties. On the other hand, 
potential models probably are also less accurate for ccc than for bbb baryons, because the ccc 
system is more relativistic and spin-dependent contributions may start playing a signihcant 
role. Although of small importance in heavy quark systems for being suppressed as Mg^, 
the spin-spin interaction derived from the OGE may help to improve the description of 
the nonperturbative lattice QGD results. One may therefore add a spin-spin term to the 
interacting potential, having the hnal form. 




i<j 


Ti 
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i<j 


Ml 


rr, 


2 ' ^3 


(3) 


We have performed a ht to the ccc data using the same Goulomb and linear conhning 
strengths used for the bbb system, A = 0.1875 and B = 0.1374 GeV^. When the spin-spin 
term is considered, one can reduce the value of the charm quark mass preserving a nice 
agreement with data. We show in Fig. [3] the results for standard values of quark-potential 


models, rric = 1.6 GeV and tq = 0.28 fm [25|]. As one can see there is a good agreement. 
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FIG. 3: ccc excited state spectra, solid lines, for the potential of Eq. ([3]) including a spin-spin 
interaction. The boxes stand for the results of Ref. I™. 



sustaining the result of a short-range Coulomb potential supplemented by a linear confining 
interaction. The short-ranged spin-spin interaction helps to the fine tuning of the negative 
parity excited states (where it is negligible) with respect to the positive parity excited states 
(where it is attractive) due to the existence of good diquarks with spin 0 in the wave function. 

One may hnally wonder if the Coulomb and linear conhning strengths of the 3Q poten¬ 
tial that we have found reproducing the nonperturbative lattice QCD results of Refs. 1,0 
are unique. To answer this question we have considered the most relevant states, the 
hrst, E2{3/2~^), and second, E^{3/2~^), radial excitations and the hrst orbital excitation, 
£^i(l/2“)^. We have represented in Fig. IHthe sets of parameters that are able to simulta¬ 
neously reproduce the positive and negative parity excited states for the bbb system with 
a reasonable bottom quark mass, mj, = 4.655 GeV. As one can see there is a small range 
of parameters reproducing the general pattern of the bbb spectrum. The strengths of the 
Coulomb and linear conhning potentials we have determined giving rise to the results shown 
in Fig. [2] (A = 0.1875 and B = 0.1374 GeV^), are in the middle of the aforementioned area 
where the most important states are properly reproduced. In Fig. [5] we have represented 
the parameters describing the bbb and the ccc systems with a simple Goulomb plus a linear 
conhning potential, using a reasonable charm quark mass, rric = 1.6 GeV. There are mainly 
two noticeable aspects. First of all, the smallness of the area reproducing the ccc data. Even 
more, such area is in the limits of the uncertainties of the nonperturbative lattice QGD data 
for the ccc system. The second relevant aspect is that there is no overlapping between the 
ccc and bbb range of parameters. As has been already noted, these results could indicate 


^ Having in mind the level degeneracies existing in the absence of rotational symmetry breaking, these 
states represent the majority of the spectra Q. 
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FIG. 4: Coulomb and linear confining strengths reproducing selected states of the bbb system. 
The light blue area stands for sets of parameters reproducing the E2{3/2~^) energy, the green 
area reproduces the £’ 2 ( 3 / 2 +) and £ 3 ( 8 / 2 +) energies, the read surface reproduces the £ 2 ( 8 / 2 +), 
£ 3 ( 3 / 2 +), and £i(l/2“) energies. The dark blue point corresponds to our best ht. 


^ccc and 



FIG. 5: The light blue area stands for Coulomb and linear confining strengths reproducing the 
£ 2 ( 3 / 2 +), £ 3 ( 3 / 2 +), and £i(l/2“) energies of the bbb system. The green surface represents the 
same for the ccc system. The dark blue point corresponds to our best ht. 
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FIG. 6: The light blue area stands for Coulomb and linear confining strengths reproducing the 
E2[2>/2'^), £ 3 ( 3 / 2 +), and £ 1 ( 1 / 2 “) energies of the bbb system. The green surface represents the 
same for the ccc system considering the spin-spin interaction of Eq. ([3|). The dark blue point 
corresponds to our best fit. 


that the spectrum of the ccc system hardly accommodates to a simple linear confining plus 
Coulomb potential, needing of higher order terms in the interaction. If we repeat the same 
simulation for the ccc system adding the spin-spin force of Eq. ([3]), we get the results shown 
in Fig. ini where one can see the overlapping between the bbb and ccc sets of parameters that 
coincides exactly in the Coulomb and linear strengths we have determined. 

Let us finally comment on the possibility of applying the approach outlined in this work to 
light flavor sectors, as for example the excited spectrum of the baryon that has also been 
studied by means of lattice QCD js^. Unfortunately the role played by other contributions 
beyond the conhnement and the Coulomb potentials, as illustrated in Ref. j3l| , leads to the 
conclusion that only baryons made of heavy flavors are clean enough to properly disentangle 
conhnement and Coulomb contributions. 


IV. SUMMARY 

To summarize, the spectra of baryons containing three identical heavy quarks, b or c, have 
been recently calculated in nonperturbative lattice QCD. The energy of the ground state 
and the lowest positive and negative parity excited states has been determined with high 
precision. These achievements constitute a unique opportunity to test phenomenological 
potential models in the regime in which they are expected to work best. We have analyzed 
these results by means of a Cornell-like potential using different numerical techniques for 
the three-body problem. For the case of the bbb system a good agreement is obtained by 
means of a simple Coulomb plus linear conhning potential. For the case of the ccc system 
the additional contribution of the spin-spin interaction is needed to have a reliable mass 
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for the charm quark. The spin-spin interaction comes suppressed by but it helps to 

correctly allocate the negative parity excitations with respect to the radial excitations of the 
3/2+ ground state. As in the case of the heavy meson spectra, a larger value of the Coulomb 
strength than predicted by SU(3) lattice QCD is needed. The phenomenological strengths 
of the Coulomb potential reproducing the heavy meson and the triply-heavy baryon spectra 
satisfy A/a < 1/2, slightly different from the 1/2 rule as the one-gluon exchange result. The 
strength of the linear conhning interaction has to be also larger than derived from SU(3) 
lattice QCD. Our results support a coherent description of the bbb and ccc heavy-baryon 
spectra with the same Coulomb and conhning strengths. 

Let us conclude by emphasizing that the recent improvements in lattice QCD calculations 
of the heavy-baryon spectra js^ may beneht of constituent quark model predictions, and, 
on the other hand, potential model approaches do also require guidance of lattice data 
for constraining not only the parameter space but also the possible functional forms to be 
explored. 
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